The IZq Approach to Epidemic-non-Epidemic Phases Revisited 

O.E. Aiello and M.A.A. da Silva 

Departamento de Fisica e Quimica da FCFRP, 
Universidade de Sao Paulo, 14040-903 Ribeirdo Preto, SP, Brazil 
(Dated: February 13, 2003; Received text; Revised text; Accepted text; Published text) 

Abstract 

In this work, we revisit the basic reproduction rate TZq definition for analysis of epidemic-non- 
epidemic phases describing the dynamics of the discrete stochastic version of the epidemic SIR 
model based on the Master Equation formalism. One shows that it is a very precise and efficient 
way to determine the epidemic threshold; using its most primitive concept, we can find exact 
results. 
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Introduction - The basic reproduction rate, TZo, is the most fundamental parameter used 
by epidemiologists |l]. It has raised interest of physiscists because of analogies between 
epidemic and percolation systems S, 0, 0] , and its generality, that permits, for example, 
analyze spreading of viruses in a structured scale-free network j^. The definition of TZq is 
"the average number of secondary cases caused by an infectious individual in a completely 
susceptible population" This simple idea had a profound effect on epidemic theory. It 
may be a global insight that cuts through the details of the transmission process, because 
it originated from consideration of deterministic models of homogeneous population with 
random mixing . However, putting more details in the model one can extend its definition 
to heterogeneous mixing. In this way, naturally the efforts were in analytical calculations of 
TZq for continuous deterministic or stochastic models 0, 0, . The stochastic framework, 
alt 

fl 

machine technology has spread the use of computationally intensive methods to solve a 



l oug h more realistic in principle, it is more complex to analyze because of detail required 
Commonly, simulations helped to confirm theoretical assumptions. Improved 



Hat diversitv of epidemics, and so simulation technics, as Dynamical Monte Carlo (DMC) 
15L lid Il7| . are becoming more popular in this subject. Taken on the advantage of 
the DMC method, we can calculate TZq straightforward. Thus, we adopted this approach 
to characterize the separation between the epidemic and non-epidemic phases, because its 
efficiency and simplicity that enables a facile finding of exact results. 

Throughout this Letter, we shall consider the classical SIR (Susceptible, Infected, Re- 
moved) epidemic model, originally based on the chemical "mass action" principle (see 3| 
and references therein) , to illustrate the discrete stochastic approaches. Our model encloses 
the deterministic one as a particular case. Based on the Master equation formalism we 
generalize the SIR model and TZq for discrete stochastic systems. Also, we describe the 
local epidemic model with an exact result to TZq using its primitive concept. Finally, we 
show a phase diagram of the local — SIR model, with deterministic-like behavior, where no 
homogeneous mixing is considered, and dynamical Monte Carlo simulation results. 

Generalized SIR model and Ro - Stochastic process approaches could simulate non- 
equilibrium systems, even the deterministic ones, introducing random variables to describe 
them in a microscopic scale. The macroscopic behavior of some system is resulting from 
averages of its microscopic properties. One can describes the evolution of the distribution 
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of probabilities, for markovian processes, with the Master Equation: 

dm) 



dt 

3 3 

where is the probabihty to find the system at the state i at the time t, and Wi^j is the 
transition probabihty per unity of time. Considering Tij the probabihty of transition from i 
to j, we may write Wi^j = ^ where Xj is a characteristic time constant (lifetime) of 

the state i. 

We now start by choosing a convenient physical extensive microscopic quantity Ai, which 
depends only of the system's state i. Since the time must change for every successful event, 
we will consider only counting events related quantities. To SIR epidemic systems the 
number of infected individuals, for example, is an adequate quantity because it represents 
the balance between the number of infection and removal events. The mean value for a given 
quantity at the time t is 

Ait) = {A) = J2m)A. (2) 

i 

This equation represents a continuous physical macroscopic quantity A{t). We can differ- 
entiate both sides of the equation above, with respect to t. After that, using ((H), and by 
defining AAij = A^ — Aj, we get 

dA{t) 



dt 



5^5^u;,^,P,AA,. (3) 
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Consider now the nearest-neighbor states j of a given state i; if we measure the "distance" 
between the states, say by the quantity |Ay4jj|, such that the non-null minimum value is 
|A74jj| = a, we may approach the equation Q by: 

ii3) 

where the symbol {ij) denotes a nearest-neighbor pair of states, and 5ij = AAij/\AAij\. 
Now we consider another physical quantity A'^ that represents a source for the quantity A. 
Thus, we can rewrite (jH) as: 

dA{t) 



dt 

3 3 



where rj =< Wj^i >i are the averaged transition probabilities per unity of time over the 
ensemble of the nearest-neighbor states i of j at some time t, i.e., the mesoscopic rates. 



Here, the word ensemble means a set of configurations accessible at a finite (small) time 

n 

around a time t; in this sense we are using a time dependent ergodicity idea and so 
generally the systems are non ergodic in nonequilibrium states. The superscripts " + " 
and " — " mean respectively the contributions to increasing and to decreasing the quantity 
A{t) |2fl]. 

Based on (jS)), we formulated the GSIR model through the following set of stochastic 
differential equations and inter-classes rates: 

f = T.''HsP3R3-J2'hPjS,, (6) 

j j 

j j 

= Y^^WP^h -Y.^'RsP3Rr (8) 

j j 

The mesoscopic rates are , and f'^jg-, for each state j, respectively, from 5* — * /, 
I R and R —>■ S. To satisfy the SIR condition the set {r^'j^s} "^^vll. Note that we 
meant that, for example, if A = J, then A'' = S in the equation (0). The conservation law 
with the total number of individuals = S(t) + I(t) + R(t) is satisfied. One may obtain 
the reproduction rate, TZq, directly from the equations (jUl— THj) with the epidemic condition 
^ > 0; where the equality is the threshold and it is set to to = 0, the initial time. Thus, 

we can do J2j''^siPj^j ~ J2j''^'iRPj^j ^ 0' what implies that J2j ^siPj^j ^ J2j''^mPj^j- One 
can thus write the reproduction rate as 

for stochastic processes. As the condition to TZq to the epidemic threshold must be valid to 
the ensemble average of initial states jo that gives the same initial condition 5*0, Iq and -Ro; 
we may define 

< ^IR >0 -'0 

the average number of the secondary cases produced by Jo infected initially. Note that if 
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we do < rsi >o= ^Iq and < rm >q= a we recover the deterministic case 
observe that if some initial configuration is fixed, one does not need the averages in (jlOj] . 
but only obtain the rates, rsi and r//j; so, of course, generally, the TZq depends on the initial 
configuration choice. In many practical situations this is important because it will determine 



an epidemic or not. We can easily adapt the result above to other models, as the 5*/ 5* model 
Q], for example, to analysis of the epidemic threshold. 

Local epidemic model - Generically, the temporal and spatial evolution characterize any 
epidemics, where in each part of the system the density of the elements can vary with the 
time. One can analyze this process through a two-dimensional lattice, in which each site, 
representing an individual of the population, receives own attributes as susceptibility and 
interactivity referring each site with the others. We will analyze a model with local contact 
only. The elements are all fixed, i.e., no populational mobility is considered. The main 
reason to study a such particular model is that we have more fluctuations, and so it is a 
good test to the efficiency of our approach. 

The probability of individuals contracts the illness, in transmitting a disease by contact, 
depends on the status (susceptible, infected or removed) in which they meet its neighbors; 
its chance of getting sick will depend on the number of sick neighbors. Thus, considering an 
element possesses n infective neighbors, and an infection chance, po, due to each neighbor, 
the probability of its change in a sick element (through n effective contacts) will be jl8l |: 

w,, = A[i~ii-pon (11) 

Therefore, (1 — po)" is the probability of no infection of a susceptible (individual) if it has n 
infected neighbors, thus 1 — (1 — Po)" is the probability of infection of a susceptible if it has 
n infected neighbors. The A parameter gives the w^^- as inverse of time units. When n = 0, 
that is, when no neighbor is contaminated, the probability of contamination due to contact is 
zero, so w^j increases when the number of effective contacts, n, increases. A global removal 
rate determines the infectious period, and an infected individual turn immune (removed) 
stochastically. So the infectious period for each individual fluctuates over an average number 



given by the inverse of the removal rate, like in the mean field approach 2l|. In this sense 
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Also has a difference 



our definition of TZq is more general than that gave in the reference 
that it is considered instantaneously instead during the infectious period, so it follows close 
the classical definition. However, we considered a range of initial values to the number of 
infected individuals, and did an analysis with an initial random distribution of immunes. 
Note that for the considered model in a square lattice, when Jq = 1, we have the trivial 
exact result 

TZo = (12) 

WIR 



where nmax is the maximum number of contacts. One can see that the exact result to a 
more thorough model, including the homogeneous mixing (mean field), is straightforward. 
The reference jisj shows this result as an analytical approximation to TZq. 

Results and final remarks - For practical purposes one distributes the individuals on a 
square lattice of N = M x M sites. All the individuals at the lattice boundary have their 
statuses fixed at susceptible status. One considered only two lattice sizes, M = 10 and 
200, because increasing M reaches smooth curves near to the M = 200 rapidly. We did 
Vjj^ = q = 1, constant independent of the configuration, and = < wsi >j, averaged over 
all individual probabilities that appear in any random configuration j; with wsi modeled 
with a purely local interaction with A = 1, so wsi = 1 — (1 — Po)"- The initial condition 
for the number of infectives, Iq, to the system, for M = 10, was varying from 1 to 99, and, 
for M = 200, from 1 up to 38000 randomly distributed on the lattice. One occupies the 
remaining sites hj Sq = N — Iq susceptibles, so -Rq = for both cases. We did the variable 
n as an integer ranging from to 8, since the first and second nearest infected neighbors 
are indistinguishably considered for each susceptible. 

Considering that the sum of the "microscopic infiuences" creates a rate rsi, we calculated 
by a sample average the initial mean rate < rsi >o- Of course in the practice of the 
simulation we drew randomly only a few configurations to estimate the averages and we 
used the equation ()10|) to elaborate the phase diagrams for the SIR model showed in the 
figure 1. We got the values for < rsi >o using 100 and 4 x 10^ configurations for M = 200 
and 10 respectively. Observe that the maximum TZq happens when the probability po reaches 
its largest value in a great population of susceptible. We found for M = 10, the epidemic 
threshold in the interval TZq = 1.03 — 1.05. For M = 200, the interval reduced to a very thin 
line defined at TZq = 1, as expected to infinite systems. Note that the TZq contour lines for 
M = 10 are noisy because of the small size of the lattice. For cases when TZq > 1, happen 
epidemic bursts in average, and when IZq < 1 does not. We did, also, experiments with a 
non-zero initial immune individuals number, Rq, and its qualitative effect on TZq is the same 
as for Rq = 0, since we have a random initial distribution of immunes. Quantitatively, for 
the initial infectives number fixed, we need to increase the probability per unit of time, wsi, 
to have epidemics; as expected, has a critical value to Rq that no epidemic occurs, the so 
known herd immunity effect [l7|. For the considered parameters of our model, we are at the 
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epidemic threshold when 

< rsi >l= ^. (13) 

After several fittings we found the reasonable expression to the initial average rate: < 
^si >o= ^ (1 ~Po)^], with the real numbers a < 1 and < (3 < 8. The critical value for 
Ro is so: Ro = N — Io[l + 1 / < rsi >q]- Note that if the condition Jq > aSo is satisfied we 
have no epidemics for any value of Pq. At first sight it might be a strange result, however, 
if the initial state has most infected individuals, the removal chance of some is high; so it 
decreases faster than the number of infections itself. 

The figure 2 show some epidemic and non epidemic cases obtained by Dynamical Monte 
Carlo. This direct method was used to find the threshold line of the phase diagrams showed 
in the figure 1. The very beginning of the process is sufficient to determine TZq, but for 
completeness we showed the epidemic curves in an extended time. 

In summary, we presented a stochastic version for the parameter TZq, based on the descrip- 
tion of the SIR model, by means of the Master equation formalism. This way, the predictive 
power TZo is transported from deterministic to the stochastic one, generalizing the concept. 
In fact, a very defined result to the threshold curve was earlier found to deterministic sys- 
tems; it is interesting that we can have this, also, for stochastic systems. It is consequence 
of the definition of TZq, whose difference with that one of already cited recent work 13| did 
decrease the threshold. Fluctuations to favor epidemics in our case when TZq < 1 are smaller 
than those favor non-epidemics, i.e., the number of non-epidemic cases prevails, even for 
finite small systems; the same happens with the opposite case. Complex geometries can be 
included in this model, since have no restrictions to the model in this sense, and the system 
geometry is an important factor to change the threshold. If the system is strongly depen- 
dent on the initial conditions, the averages are not appropriated to predict an epidemic. We 
believe that this definition of TZq open the doors for new investigations and calculations of 
TZq for more realistic systems because we used a general microscopic description to get a 
parameter of macroscopic nature. 

The authors gratefully acknowledge funding support from FAPESP Grant n. 00/11635-7 
and 97/03575-0. The authors would also like to thank Drs. A. Caliri and V.J. Haas for 
many stimulating discussions and suggestions. 
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FIGURE CAPTIONS 



Figure 1. Shows a phase diagram for the local — SIR model with parameters g = 1.0, ^ = 
(0.05 — 0.95) and po — (0.05 — 0.95). The values of TZq larger than 1 in the smooth contour 
hues (M = 200) allow that the number of infected increases with the time characterizing 
an epidemic outbreak, for TZo smaller than 1 the infection fade-out. The threshold to the 
noisy contour map (M = 10) is in the interval 1.03 — 1.05 that is too small to show in the 
diagram. 

Figure 2. Epidemic curves. Show the number of Infectives evolving with the time. The 
numerical values for the model parameters are q = 1, M = 200 and /(O) = 5000. A total of 
20 experiments was done to get the averages. The figure shows two cases, some curves those 
represent epidemic outbreaks (T^-o > 1) and others in that the infection fade-outs (T^-o < 1). 
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Figure 2 
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